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ABSTRACT 

We present results of 2 mm observations of the Crab Nebula, obtained using the Goddard-IRAM 
Superconducting 2 Millimeter Observer (GISMO) bolometer camera on the TRAM 30 m telescope. 
Additional 3.3 mm observations with the MUSTANG bolometer array on the Green Bank Telescope 
are also presented. The integrated 2 mm flux density of the Crab Nebula provides no evidence for 
the emergence of a second synchrotron component that has been proposed. It is consistent with the 
radio power law spectrum, extrapolated up to a break frequency of log(i^f,[GHz]) = 2.84 ± 0.29 or 
I'b = QQ^ttll GHz. The Crab Nebula is well-resolved by the - 16'.'7 beam (FWHM) of GISMO. 
Comparison to radio data at comparable spatial resolution enables us to confirm significant spatial 
variation of the spectral index between 21 cm and 2 mm. The main effect is a spectral flattening in 
the inner region of the Crab Nebula, correlated with the toroidal structure at the center of the nebula 
that is prominent in the near~IR through X-ray regime. 

Subject headings: ISM: individual (Crab Nebula) — ISM: supernova remnants — Radiation mecha- 
nisms: non-thermal 



1. INTRODUCTION 

The Crab Nebula is the prototypical example of 
a pulsar wind nebula (P WN) , plerion, or filled- 
center supernova remnant ((Davidson fc FesenI 119851 : 
IGaensler fc Sla^ [2006t IHested 120081) . As one of the 
brighter objects in the sky across the entire electro- 
magnetic spectrum from the radio to gamma rays, the 
Crab Nebula is often one of the first targets observed 
by new instruments, whether they provide new spectral 
windows, or offer improvements in sensitivity, or spa- 
tial or spectra l resol u tion (e.g. [B olton & Stanlev 194ft 
Marsden ct al. 1984; Temi m et al.l 120061: IWaller et all 
r995;.Bowver et al...l964: , Fazio et al.lll972| ). 
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At radio wavelengths the Crab Nebula's integrated 
emission is well known to exhibit a power law spec- 
trum Si, = Sqv" , with a spectral index a — —0.299 
(jBaars et al.lll977[ ). This spectral index is fairly typical 
of the synchrotron emission exhibited by PWNe, and is 
distinctly flatter than the spectral indices of more com- 
mon shell-like SNRs where the relativistic electrons are 
accelerated in the shocks of the expanding blast wave 
rather than by a pulsar wind. In the optical regime, 
the sync hrotron emission exhibits a s teeper spectral in- 
dex ('e.g. IVeron-Cettv fc Woitie^l 19931) . Thus, a spectral 
break has been in ferred to lie at wave l engths b etween 
10 and 1000 ^m (Ma rsden et all [1981 IWoltieil [1987). 
However, the sharpness of the break and the exact loca- 
tion are somewhat obscured by the difficulty in obtaining 
high precision measurements at mm wavelengths, and the 
presence of stronger th ermal emission from dust a.t mid 
IR wavelengths ( Marsd en et"alll 19841: iMezger et allll986l: 
I Strom fc Greidanusi 11993 ). The dust is associated with 
the knots and filaments of ejecta in the Crab Nebula, 
although spectral observations have not revealed char- 
acteristic features such as silicate or PAH ba nds in the 
spectrum (|Douvion et al.ll2001t l^mim et al.l l2006). 

Resolved imaging of the Crab Nebula indicates that 
there is very little variat ion in the radio spect ral index 
across the Crab Nebula ()Bietenholz et al .' 1997). There 
has been some indication of a spectral index variation 
between the radio and mm regimes, suggesting the possi- 
bility of_a_£hy^kaJjxdisti^^ second synchrotron compo - 
nent (jBandiera et al.ll2002D . although l&een et al.l ()2004f ) 
argue that the variations may in fact be due to known 
temporal variability of the SNR. 

In this paper, we further probe the synchrotron com- 
ponent using a map we obtained with the GISMO 2 mm 
camera on the IRAM 30 m telescope. The paper is 
organized as follows: in section [5] we introduce the 
various maps that were used. In the process we de- 
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scribe the GISMO instrument, the data reduction pack- 
age (CRUSH) and the data analysis. The 21 cm and 
6 cm radio maps used for comparison are introduced. In 
section [3l we study the integrated radio/sub-milhmeter 
spectrum of the SNR. We derive the location of the 
break in the spectrum. In section |4l we present and dis- 
cuss spectral index maps created by comparison of the 
GISMO data with the radio maps. We also discuss the 
correlation of features in the spectral index maps with 
features at other wavelengths (section [5]) . Finally our 
conclusions are stated in section |6l 

2. OBSERVATIONS 
2.1. GISMO Data 

The Crab Nebula was observed at the 30 m telescope 
located on Pico Veleta (near Granada, Spain) and oper- 
ated by t he Institute de Ra dioastronomie Millimetrique 
(IRAM) (jBaars et al.lll987f) . The observations were cen- 
tered at (aj2ooo,^J2ooo) = (5''37'"45^ 22°1'0"). The 
data were obtained during the first run of the Goddard- 
IRAM Superconducting 2 Millimeter Observer (GISMO) 
in November of 2007. 

The GISMO instrument is a bolometric camera devel- 
oped for operation in the 2 mm atmospheric window. 
The central frequency of the GISMO band is 150 GHz, 
with 5v/v = 0.15 (or a bandwidth of 22 GHz @ 150 
GHz). The instrument uses an 8 x 16 planar array of mul- 
tiplexed superconducting transition edge sensor (TES) 
bolometers which in corporates the Ba ckshort Under Grid 
(BUG) architecture (I Allen et al.ll2006D. The size and sen- 
sitivity of the detector array |Staguhn eFaH [20061) give 
this instrument significantly greater mapping speed at 
higher angular resolution in this wavelength than has 
previously been possible. GISMO has a pixel size of 
2 mm which provides for an angular separation of 14" 
on the sky corresponding to a sampling of 0.9A/Z? at 
2 mm wavelength. With an array of 8 x 16 pixels, the 
GISMO instrument has a gapless field of view of 2' x 4' 
which allows it to take advantage of a significant frac- 
tion of the telescope's optical capabilities. The FWHM 
of the beam at 2 mm is ~ 16'.'7, which i s ^ 15% wider 
than t he ideal diffraction limited beam. iStaguhn et al.l 
(12001 gives a more detailed description of the instru- 
ment. On this first run of the instrument, there was 
a mechanical problem that caused one quadrant of the 
pixels to malfunction. The data reduction also identified 
pixels that exhibited increased noise and these too were 
excised, leaving 62 detector pixels used for this analysis. 

The maps were obtained on November 12, 2007 us- 
ing five pairs of on-the-fly scans. The first scan in each 
pair used 33 scan lines in the azimuth, 630" long, off- 
set by 18". The second scan of each pair transposes the 
same pattern with 33 scan lines in elevation (Fig. [1]). 
Each scan lasts 512 seconds, and a pair of scans requires 
1035 seconds. The integration time per 3" pixel in the 
final map is ^7.4 sec for each pair of scans. During a 
complete raster scan each element of the bolometer ar- 
ray then samples the whole of the Crab w hich has an 
apparent size of 420" x 290" ()Trimblelll973l ). 

The data were reduced using the Com prehensive Re - 
duction Utility for SHARC-2 (CRUSH: iKwIci [20081 ) . 
CRUSH was initially developed for data reduction of 
SHARC-2 observations, but it has been enhanced to re- 
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Fig. 1. — The scan pattern for one of the five pairs of scans used 
on the Crab Nebula. An on-the-fly (OTF) scan map is executed in 
azimuth, immediately followed by another OTF map in elevation. 
The pattern of these scans, when projected into equatorial coordi- 
nates as shown here, rotates over time as the parallactic angle of 
the target changes. The parallactic angle for this pair of scans is 
~ 45°. 



duce data from other instruments, including GISMO. 
CRUSH uses an iterated sequence of statistical esti- 
mators to separate source, atmosphere and instrument 
signals. As part of this reduction, the final map is 
smoothed by a Gaussian function with 2/3 beam width 
(the CRUSH default for faint sources) to have an effective 
FWHM of ^ 20". This smoothing is done to minimize 
spurious high spatial frequency features that are related 
to the relatively coarse samp ling of the be am with the 
the GISMO detector pixels (jKovacsl [200l . The final 
map produced by CRUSH has a pixel scale of 3" (Figure 
[2]). The noise level of the map varies as a function of 
position due to the coverage of the scans, but over the 
region where the Crab Nebula is brighter than 0.5 Jy 
beam^^ the noise level is ct < 0.016 Jy beam~^. 

The total integrated flux density of the Crab Nebula is 
determined to be 244 ± 24 Jy at 2 mm, after adjustment 
for the slightly negative local background level in the 
image produced by CRUSH. This value is based on use of 
Saturn as a calibrator, which was observed after the Crab 
Nebula during a period of poorer weather (risocffz ^ 0.3 
vs. Ti50GHz ^ 0.12). 

2.2. MUSTANG Data 

Additional observations of the Crab Nebula were 
carried out with Multiplexed Squid TES A rray at 
Ninety GHz (MUSTANG; Di cker et al.l [2508l) instru- 
ment on the Green Bank Telescope (GBT). The data 
were acquired during Feb 2008, under project code 
AGBT08A_ 056. The data were calibr ated rel a tive t o 
Mars using IWeiland et all ([20Tl[) and IWrightl ([2001 . 
A map was made with a maximum likelihood imaging 
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Fig. 2.— (top) GISMO map of the Crab Nebula at 2.0 mm. The 
image is shghtly smoothed from the intrinsic 16'.' 7 (FWHM) resolu- 
tion. Across the region of the Crab Nebula, the uncertainty in the 
intensity is 0.008 < cr < 0.016 Jy beam"!, (bottom) MUSTANG 
map of the Crab Nebula at 3.3 mm. The beam is 9" FWHM. 
The gray scale bars are labeled in Jy beam~^. The axes indicate 
right ascension and declination offsets from the pulsar location; 
(aj2000,<5j200o) = (05'' 34™ 32^ +22° 00' 52"). 



pipelin e develo ped for the Atacama Cosmology Telescope 
(ACT: iFowleT e t al. 2010) that was adapted for use with 
MUSTANG (Figure E]). The MUSTANG camera has 9" 
(FWHM) spatial resolution, compared to 16'.'7 for the 
GISMO instrument. Due to the comparatively small in- 
stantaneous FOV of MUSTANG (40" x 40"), the total 
integrated flux density for the Crab Nebula is not well 



constrained. Structures on angular scales from 4" to 1.'5 
are recovered with good fidelity. 

2.3. Radio Data 

To understand the synchrotron emission from the Crab 
Nebula, we compared the GISMO 2 mm map to radio 
maps at lower frequencies. Radio maps at two different 
frequencies and two different epochs were used for com- 
parison. The first map at 1.41 GHz (21 cm) was obtained 
using all 4 c onfigurations of the VLA over t he period from 
1987-1988 ()Bietenholz fc Kronbertd IT990t) . This image 
has a cleaned beam size of 1"8 x 2'.'0. A second VLA 
map at 5 GHz (6 cm) was obtained during 2001, and has 
a cleaned beam size of 1"4 (Bictenholz ct al. 2001). The 
5 GHz image (only) has been corrected for the VLA's 
primary beam response (~ 9' FWHM at 5 GHz). 

For comparison with the GISMO image, the radio 
images were rescaled in both size and brightness. We 
have used 1987.49 as the epoch of the 1.4 GHz image, 
and 2007.87 is the epoch for the GISMO image. This 
20.38 yr time difference results in a 3.1% increase in the 
siz e the nebula acc o rding the the expansion measured 
by iBietenholz et al.l ()1991fl . This expansion corresponds 
to ~ 0.5 of the GISMO beam width at the outermost 
edges of the SNR. There is a 3.5% decrease in the radio 
brightness between the dates of the 1.41 GHz VLA and 
2 mm GISMO observations, assuming a fading rate of 
-0.167 % vr-i (lAller fc Revnoldd[l985l). M ore recent re- 
ports (|Vinvaikinll2007l: iWeiland et al.ll2011[ ) confirm that 
the fading rate is not a strong function of frequency or 
time. After applying these corrections to the radio im- 
age, the map was convolved with a Gaussian function to 
match the 20"07 resolution of the GISMO data. Equiv- 
alent adjustments were made to the 5 GHz VLA image, 
resulting in a size increase of 1.05% and a brightness de- 
crease of 1.169%. The similarity of results derived from 
both the 1.41 and 5 GHz images indicates that the cor- 
rections for expansion are sufficiently accurate at GISMO 
resolution. 

3. INTEGRATED FLUX DENSITY AND SPECTRAL INDEX 

Figure [3] shows the integrated spectrum of the Crab 
Nebula including the 2 mm GISMO fiux density, and 
other measurements from the radio to optical wave- 
lengths (30 cm - 3000 A or 1 - 10*^ GHz). All flux 
densities have been corrected to the GISMO epoch of 
2007.87. The 2 mm GISMO fiux density is consistent 
with an extrapolation of the 1-35 GHz radio spectrum 
wi th a power law inde x a = — 0.299±0.009 as determined 
bv lBaars etall (|1977D . 

At 10 /xm to 200 /xm, the spectrum exhibits a bump due 
to thermal emission from dust. Shorter wavelengths ap- 
pear to resume a power law spectrum, but with a steeper 
slope than at the radio - mm wavelengths. We charac- 
terize the spectrum by fitting a broken power law of: 

5(^) = (^('^"K'^/'^''):: ^<'^^ (1) 

characterized by the four parameters S(vi,), vi,, ai, and 
a2- The 10 - 100 /Ltm measurements are not used to 
constrain the fit because of the possible influence of the 
dust emission at these wavelengths. At the shorter wave- 
lengths we choose to constrain the fit only with the 
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Fig. 3. — The integrated radio — optical spectrum of the Crab Nebula. The 2 mm GISMO result (circled) helps confirm that the 
break in the power law spectrum at radio wavelengths must lie a wavelengths shorter than 2 mm. The line shows a broken power law 
fit to the data, with dotted regions indicating ranges where data were ignored because of thermal emission from dust (mid-IR), or the 
necessity for larger extinction corrections (optical ncar-IR). Numbers along the bottom of the plot indicate the sources of the plotted data: 
(1) Baars ct al. (1977) and references therein, (2) Macias-Pcrez ct al. (2010), (3) Weiland ct al. (2011), (4) this work, (5) Bandiera et aT] 
(feooal. fel.Mczger et al., tl98Q). (7) .Wright et al.. (1979) . (8) .Green et al.i (,2004) . (9) iTemim et al.i (i200& ). (A) iMarsden et al.i (,1984 ). (B) 
IGrasdalenI (119791) . andTC) IVeron-Cettv fc Woltiea TT9§a) . 



3.6 - 8 fj,m Spitzer IRAC measurements of iTemim et al.l 
(|2006D . because the near-IR and optical measurements 
are much more dependent on the application of ex- 
tinction corrections, and may also be affected by ad- 
ditional intrinsic steepening of the spectrum. We find 
\og[S{vb)] - 2.14±0.09 (or S{vb) = 137^11 Jy)> ^og{vb) = 
2.84 ± 0.29 (or i/f, = egstggg GHz), ai = -0.30 ± 0.01, 
and a2 — —0.50 ± 0.03. T he low fr equency spectral in- 
dex, ai, is consistent with Baar s et al . (1977), while the 
high frequency spectral index, 02, is consistent with re - 
suh (a = -0.50±0.10) obtained bv rPouvion et all (|200l 
using ISO ISOC AM data covering the central 3' x 3' of the 
SNR at 4.5 and 11.4 ^m (LWl and LW8 filters). If the 
near-IR data of Graselden were included in the fit, then 
log[S'(z^f,)] would decrease by ~ Icr, a2 would increase 
(flatten) by la, and ai would remain unchanged. Ad- 
ditional inclusion of the optical measurements would re- 
sult in very significant changes in the derived break fre- 
quency (log[S'(i'b)] ~ 3.85) and a steeper high frequency 
spectral index (a2 ~ —0.64). However this fit systemat- 
ically overshoots the IRAC measurements, especially at 
8 fim where dust and line emission should result in an 
enhancement over the synchrotron emission. Therefore 
we conclude that there is either a problem with the ex- 



tinction correction or relative calibration of the optical 
observations, or the intrinsic synchrotron spectrum has 
curvature or additional breaks such that it cannot be 
properly characterized by a single power law at mid-IR 
to optical wavelengths. 

A detailed account of the synchrotron spectrum of 
the Crab Nebula (and other pulsar wind nebulae) is 
provided by the diffusive synchro t ron ra diation (DSR) 
model of [Fleishman fc Bietenhoig (|2007[ ). In the DSR 
model, the magnetic field is tangled on a range of spa- 
tial scales characterized by a power law index (n) and 
a maximum scale length (L). The predicted emission 
spectrum will contain a break at a frequency that is 
a function of n and L, the magnetic field strength 
(B), and the minimum energy (7 = E/rac^) of the 
power law distribution of relativistic particles, specif- 
ically Vb (X j^(i-»)/(i-Hn)^2/(i- f-n)^2_ Adopting nomi- 
nal va lues of these parameters [Fleishman fc Bi etenhold 
(|2007D find the break frequency to be at Vb « 3200 GHz. 
A shift in the break frequency may be accommodated 
by adjusting any or all of the parameters. With n being 
fairly strictly constrained to values of 1.54 < n < 1.6 by 
the spectral index of the low frequency (radio) emission 
(n = 1 — 2ai), the value of Vb will be most sensitive to 
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changes in 7 and least sensitive to changes in L. 

The GISMO r neasurement and o ur analysis support 
the conclusion of iGreen et alj (|2004D that the integrated 
spectrum of the Crab Nebula continues as a power law 
from the radio regi me down to w avelengths < 850 /im. 
While the BandicraeTaD (|2002i ) 1.3 /im measurement 
lies slightly above the power law, the integrated spectrum 
including newer data does not exhibit a significant excess 
component at millimeter wavelen gths. The 170 /im ISO 
measurement ijGreen "eralll200l seems to confirm the 
location of the break. Future observations in the 200 - 
800 /im range should help to better refine the location 
and sharpness of the break. 

4. SPECTRAL INDEX MAPS 

Maps at 20" resolution of the Crab Nebula's spec- 
tral index between 2 mm and radio wavelengths are 
shown in Figure 01 These maps show significant struc- 
ture in the spectral in dex, and strongly rese mble those of 
iBandiera et all pOOl (1.3 mm - 21 cm) and Gree n et al.l 
(|2004D (85(^7™ - 21 cm). There is a slight offset be- 
tween the mean spectral index of the two maps because 
the radio data themselves are not an exact match to the 
a — —0.299 spectral index. A multiplicative scale error 
in the absolute calibration of the radio (or GISMO) data 
will affect the mean spectral index, but not any relative 
variations in spectral index. Therefore, the contours cho- 
sen for Fig. Ulare slightly offset to emphasize the strong 
similarity in structure between the two spectral index 
maps. 

The random uncertainties (noise) in the intensity maps 
propagate into uncertainties in the derived spectral in- 
dex. For a = ln[S'(i^i)/S'(i^2)]/ln(i^i/i^2), we have 
= {[a5(.,)/5(z.i)]2-f [a5(.,)/5(^2)]n°-Vln(^i/^2). 
For spectral indices involving the GISMO (150 GHz) 
and VLA (1.4 or 5 GHz) data, the random uncertain- 
ties are dominated by the GISMO data. Except at 
the faint edges of the SNR, the random uncertainties 
in the GISMO data are < 3% (See section 2.1). Thus 
< {as(i5o)/5(150)}/ ln(150/5) « 0.01. 

Although the random uncertainties spectral index 
maps are relatively modest, there are at least four po- 
tential sources of systematic errors which may affect the 
apparent structure in these spectral index maps. The 
first is the significant difference in the resolution of the 
original maps. The radio maps were convolved to the 
resolution of the GISMO image fi i2.3p . However, if this 
process did not accurately reproduce the GISMO beam, 
then there would be artificial changes in the spectral in- 
dex in regions where there are steep brightness gradients. 
To check for errors introduced by a resolution mismatch, 
the convolved resolution of the radio map was changed 
by a factor of 50%. This produced minimal changes the 
spectral index map, indicating that the observed spectral 
index variations are not a result of mismatched resolu- 
tion. 

A second potential source of error is improper back- 
ground subtraction. The background levels were checked 
by examining the correlation of the radio and 2 mm in- 
tensities at intermediate brightness levels (100 < /radio < 
400 MJy sr~^). The correlations between the radio and 
GISMO images are very linear and extrapolate to zero in- 
tercepts indicating 2 mm background errors of -|-2.2 MJy 
sr~^ with respect to the 21 cm data, and —2.2 MJy sr~^ 
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Fig. 4. — Radio - 2 mm spectral index maps, (top) Map of the 
spectral index a, for the 21 cm to 2 mm range. The contour levels 
are drawn at a = —0.30, -0.28, -0.26, and -0.24 (white contour), 
(bottom) Map of the spectral index a, for the 6 cm to 2 mm range. 
The contour levels are drawn at o = —0.325, -0.305, -0.285, and 
-0.265 (white contour). The small offset in contour levels is done 
to highlight the similarity in structure of these two images, rather 
than the small shift in the mean spectral index found in these 
comparisons to two independent radio images. The axes indicate 
right ascension and declination offsets from the pulsar location: 
("J2000,'5j2000) = (05*^ 34™ 32^-H22° 00' 52"). 



with respect to the 6 cm data. Comparison with Figure 
[2] indicates that except in the outermost portions of the 
SNR, this background error is < 10% of the observed sur- 
face brightness, and thus any induced errors in the spec- 
tral index should be < 0.03, (much less in the brighter 
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central regions). Furthermore, as the background errors 
have opposite signs with respect to the 21 and 6 cm im- 
ages, yet the spectral index maps are very similar, we 
conclude that improper background subtraction cannot 
account for the observed structure in the central regions 
of the spectral index maps. However, at the edges of the 
Crab Nebula, the emission of the SNR falls to levels at or 
below the uncertainty in the background. Therefore, we 
cannot reliably determine changes in the spectral index 
at the edges of the SNR and do not draw conclusions 
about the spectral index in the outermost ~ 1' of the 
maps shown in Figure |4l 

A third possible source of error is that either the radio 
or the 2 mm images may be missing flux at particular 
spatial scales. This may clearly affect the radio data due 
to the limited u — v coverage provide by an interferom- 
eter such as the VLA. It may also affect the GISMO 
data due to the periodic nature of the scan pattern and 
the processing needed to separate and remove the tem- 
poral variation of the atmospheric emission. The radio 
maps were created using Maximum Entropy Deconvo- 
lution provided with a low- resolution model to account 
for flux at missing spatial scales. The GISMO maps were 
reduced using the CRUSH package which iteratively re- 
moves celestial, foreground, and instrument components 
as measured with the scan pattern one-by-one until all 
that remains is the source and white noise. Figure [5] 
shows the normalized spatial power spectra {P{k) as a 
function of spatial frequency k) of the 2 mm GISMO 
map, the 3.3 mm map, and the 21 cm radio map. For 
each data set, the power spectrum is calculated on the 
images at their original spatial resolution. The 2 mm and 
21 cm power spectra are well correlated at low spatial fre- 
quencies (fc < 0.02, or spatial scales > 50"). Because of 
the larger PSF, the power seen by GISMO (and MUS- 
TANG) becomes increasingly attenuated with respect to 
that of measured by the VLA at smaller scales. At the 
smallest angular scales, the MUSTANG power spectrum 
flattens due to white noise (instrumental in nature, and 
thus unaffected by the beam). This comparison indicates 
that missing flux at certain spatial scales is not a likely 
source of the observed spectral index variations. 

The fourth potential source of error may be caused 
by the difference in the epochs of the two maps. 
In addition to the overall evolution in the Crab 
Nebula's size and brightness, there have been ob- 
servations of very rapid changes at small spatial 
scales at op t ical, X-ray, and radio wavelength s (e.g . 
Hester et al] Il995t iGreiveldinger fc Aschenbac"bl 119991: 
Bietenholz et al.ll2004[) . However, the similaritv between 
the 21 cm and 6 cm spectral index maps provides a strong 
indicator that the temporal variations are not significant. 
For temporal changes to be responsible for the observed 
spectral index variations, they would need to be negli- 
gible between 1987 and 2001 (the epochs of the radio 
maps), and then rapidly become significant in the inter- 
val from 2001 to 2008. 

We also attempted to calculate the 2-3.3 mm and 
1.3-2 mm spectral index maps using the GISMO data 
and either the 3.3 mm MUSTA NG data or the 1.3 mm 
MAMBO ()Bandiera et al.ll2002[ ) observations. However, 
because of the close proximity of the wavelengths, small 
errors in the background or large scale structure lead 
to magnified effects in the derived spectral index, which 
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Fig. 5. — The normalized power spectra of the 2 mm GISMO 
map, the 3.3 mm MUSTANG map, and the 21 cm (1.4 GHz) radio 
map, each at its full spatial resolution. At lower spatial frequencies 
(fe < 0.02 arcsec"^) the power spectra are in fair agreement. The 
larger beams of the GISMO and MUSTANG data cause some at- 
tenuation with respect to the VLA data at the higher frequencies. 
White noise flattens the spectra on scales smaller than the beam 
(e.g. MUSTANG data). 

cause the resulting spectral index map to be unreliable. 
The mm data sets alone are insufficient to provide accu- 
rate spectral index maps over such small ranges in wave- 
length. 

5. DISCUSSION 

Having ruled out various sources of systematic error 
above, we conclude that the observed spectral index vari- 
ations shown in Figure |4] are real properties of the SNR, 
and are not caused by systematic effects inherent in the 
data or the analysis procedures. The reality of the spec- 
tral index variations is streng thened by a ve r y stro ng 
correlation with those found bylBandiera et aTl ( 20021 ) at 
1.3 mm and lGreen et al.l (|2004[ ) at 850 /im using different 
instruments. Additionally, in Figure [6] we compare the 
spectral index with the synchrotron emission as depicted 
at IR and X-ray wavelengths. The regions of flatter spec- 
tral index are well correlated with the toroidal structure 
of emission that surrounds the pulsar. This structure is 
present at IR wavelengths, but is most distinct at X-ray 
w avelengths. 

iBandiera et al.l ()2002| ) have argued that the spectral 
index variations, along with the integrated 1.3 mm flux 
density which appeared to lie above the extrapolation 
of the radio synchrotron spectrum, were indications of 
a distinct synchrotron component generated by a sep- 
arate population of r elativistic electrons . With mea- 
surements at 850 /im, iGreen et al.l (|2004f ) found similar 
spectral index variations, but with better signal to noise 
and correspondingly less apparent small-scale structure. 
However, their integrated flux density is a good match 
to the extrapolation of the radio spectrum. They con- 
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Fig. 6. — Comparison of the spectral index variations witli synchrotron emission at IR and X-ray wavelengths, (top left) Spitzer IRAC 
4.5 /im image of the Crab Nebula. Synchrotron emission is dominant at this wavelength, (bottom left) Chandra ACIS X-ray image of the 
central portion of the Crab Nebula. The toroidal structure around the pulsar is very clear at this wavelength, (top and bottom right) 
Spectral index contours a, for the 21.0 cm— 2.0 mm range (from Fig. (JJ overlaid on the left hand images. The areas of flattest spectral index 
are correlated with the brightest part of the toroidal structure surrounding the Crab Pulsar. This is especially clear at X-ray wavelengths. 
The axes indicate right ascension and declination offsets from the pulsar location; (aj2000i '5j200o) = (OS'' 34™ 32'',-|-22° 00' 52"). 



elude that a second synchrotron component is not re- 
quired, and suggest that the known temporal variations 
in the region near the pulsar may be responsible for the 
flatter s pectral index i n the region. From analysis of 
WMAP (Page ct al. 2007), Archcops (De sert et all2008l : 
iMacfas-Perez et al..,2007.'). and publi shed integrated flux 
densities, acfas-Perez et al.l (|2010() also conclude that 
the data are not significantly better fit by the inclusion 
of an additional synchrotron or cold dust component to 
enhance emission at mm wavelengths. Our integrated 2 
mm flux de nsity measured by GISMO supp orts the con- 
clusio ns of iMacias-Perez et al] ()20I0f ) and iGreen et al.1 
(|2004f ). and weaken the argument for a second compo- 
nent. However, because we find the same spectral index 
variations when comparing to radio images for two dif- 
ferent epochs, we do not believe that temporal variation 



between the times of the radio and mm observations can 
be responsible for the spectral index variations. 

The presence of small-scale spectral index variations at 
wavelengths A < 2 mm implies the integrated spectrum 
of the Crab Nebula is probably better characterized by a 
more gradual change in slope rather than a sharp break 
from one spectral index to another (Eq. 1). However, 
the present far-IR to mm data are insufficient to reveal 
this distinction in the integrated emission. The spectral 
index variations of the synchrotron emission also impact 
study of the dust emission, which is strongest between 
10 and 100 fim. Accurate spatial and spectral informa- 
tion for the dust requires subtraction of the synchrotron 
component. The synchrotron emission cannot be simply 
extrapolated from longer or shorter wavelengths using 
the global spectral index. If local variations are not in- 
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eluded, residual synchrotron emission will mistakenly be 
attributed to dust, which will adversely impact assess- 
ments of the dust temperature, mass, composition, and 
spatial distribution. 



6. CONCLUSION 

We present 2 mm observations of the Crab Nebula us- 
ing the Goddard IRAM Superconducting 2 mm Observer 
(GISMO) at the IRAM 30 m telescope. As these are 
among the first GISMO observations, they provide a test 
of the instrument, the observing strategy, and the data 
reduction with respect to extended sources. 

The GISMO map for the Crab Nebula gives a total 
integrated flux density of 244 ± 24 Jy which is in accor- 
dance with the known radio/sub millimeter spectrum of 
the Crab Nebula. The radio end of the spectrum is de- 
scribed well using a power law with ai = —0.30 ± 0.01 
which matches the well-established value of iBaars et al.l 
(|1977f ). A simultaneous power law fit to synchrotron 
emission at 3 - 8 /im allows us to locate the break in the 
spectrum to be at = 695^335 GHz. However in detail, 
the break may not be sharp, and the frequency may vary 
at different locations within the nebula. We observe no 
excess of flux which would indicate the prese nce of a new 
synchrotron component as proposed in Bandiera et al.l 
(|2002). 

Spectral index maps were created to compare the 
GISMO data with the radio regime. Two VLA maps at 
1.4 GHz and 5 GHz were used. Comparisons with these 
maps produce similar results and both seem to show a 
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